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It is well established that for non-interacting 
electrons, increasing disorder drives a metal into a 
gapless localized Anderson insulator [1 J. While in 
three dimensions a threshold in disorder must be 
crossed for the transition, in two dimensions and 
lower, the smallest amount of disorder destabi- 
lizes the metal. The nature of the metal-insulator 
transition in an interacting system remains unre- 
solved [iSHS]. Here we explore the effect of disor- 
der on a strongly correlated Mott insulator with- 
out changing the carrier concentration. Angle 
resolved photoemission spectroscopy (ARPES) 
measurements on copper intercalated single crys- 
tals of the layered dichalcogenide lT-TaS2 reveal 
the presence of new delocalized states within the 
Mott gap. This is the first experimental realiza- 
tion of a novel disorder-induced metal that was 
theoretically predicted [9j to exist between the 
Mott insulator and Anderson insulator. 
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FIG. 1: Schematic figure showing the effect of disordering a 
Mott insulator (without adding carriers). V is the disorder 
scale relative to the bandwidth. The charge gap A gap in the 
Mott insulator vanishes at the quantum phase transition to 
a metal. The novel metal has a finite dc conductivity adc 
and other unusual spectral properties elucidated by ARPES 
experiments reported here. At higher disorder we expect a 
second quantum phase transition from the metal to a gapless 
localized insulator. If the Mott insulator has magnetic order, 
characterized by the order parameter m^, it need not be tied 
to the metal insulator transitions. 

The transition-metal dichalcogenide lT-TaS2 is a 
quasi-2D system with strong electron-electron and 
electron-phonon coupling with a very rich phase diagram. 
Upon decreasing the temperature, it undergoes several 
CDW phase transitions: from a high temperature metal- 



lic phase to an incommensurate CDW (ICCDW), to a 
nearly commensurate CDW (NCDW), and finally to a 
commensurate CDW (CCDW). It even becomes super- 
conducting when subjected to external pressure or chem- 
ical disorder fiO\ [TlJ. 

Despite the seemingly complicated electronic proper- 
ties, lT-TaS2 has a simple crystal structure composed of 
weakly coupled layers, each layer containing a single sheet 
of tantalum atoms, sandwiched in between two sheets of 
sulfur atoms. The Ta atoms within each layer form a 2D 
hexagonal lattice. The basic CDW instability is formed 
within the tantalum layers by the arrangement of 13 Ta 
atoms into a "star-of-David" shaped cluster, where 12 
Ta atoms move slightly inward towards the 13th central 
Ta atom (Fig.|2|a)). As temperature is lowered, the size 
of the CDW domains becomes larger until all the do- 
mains interlock into a single coherent CDW modulation 
extending throughout the layer in the CCDW phase be- 
low TccDW=^^^^' The new "star-of-David" unit cells 
form a hexagonal lattice and the corresponding reduced 
Brillouin zone is depicted by the small green hexagon in 
Fig-Htb). 

The transition into the CCDW phase is accompanied 
by a metal-insulator (MI) transition into a Mott phase 
[ij |2] . The insulating behavior can be clearly seen from 
resistivity measurements shown in Fig.[2|g). It has been 
shown [3] that out of the 13 Ta 5d electrons within the 
"star-of-david" cluster 12 occupy states below the Fermi 
energy, while the 13th one occupies a half- filled band 
crossing the Fermi energy. This system is widely regarded 
as a realization of a half-filled 2D Hubbard model on a 
triangular lattice where strong electron-electron interac- 
tion opens a Mott gap [2]. 

In this paper, we show that even in the presence of 
disorder due to Cu intercalation, the CCDW phase exists 
below TccDW and that the low energy physics is well 
described by a disordered half-filled Hubbard model. We 
therefore have a candidate system for investigating the 
effect of disorder on a Mott state. 

Fig[2|c-d) show the Transmission Electron Microscopy 
(TEM) images for disordered lT-TaS2 above and below 
T^CCDW' The locking of the "star-of-David" clusters re- 
sults in the appearance of new points in the diffraction 
image. These points form a small hexagon (green in the 
figure) rotated by about 13^ relative to the large hexagon 
which represents the basic trigonal unit cell (red in the 
figure) . The results are identical to that of a clean sample 
(not shown here). 

The transition to the CCDW phase can also be seen 
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FIG. 2: Mott phase in lT-TaS2: (a) visualization of the CCDW phase within a single tantalum sheet and the arrangement 
of the star-of-David clusters into a triangular lattice, (b) The 2D Brillouin zone of the original triangular Ta sheet (red), and 
the CCDW phase reconstructed reciprocal unit cells (green), (c-d) TEM diffraction images from a disordered sample in the 
CCDW phase (T=123K), and the NNCDW phase (T = 223K). (e-f) ARPES data from a disordered sample showing the band 
dispersion along the F-M direction, in the NCCDW and CCDW phases, respectively. The dashed lines are guide to the eye 
showing the band dispersion, (g) Resistance as a function of temperature on for a clean and a disordered sample, the sample 
were cooled from 300K down to 2K and back to 300K. 



directly using APRES. In Fig ^e-f) we show ARPES 
(experimental details are given in the Methods section) 
spectra taken above and below the transition. Above the 
transition we observe a parabolic band with a bandwidth 
of about ^ leV, originating from the Ta 5d orbital [15], 
dispersing upwards from M towards the F point, and 
crossing the Fermi energy at 25 % of the F — M distance 
(dashed line Fig|2|e)). At 25K this band is split into 
three sub-bands due to the CCDW reconstruction (Fig 
[2|f)). The same reconstruction happens in the clean case 
also (Fig. [sj^a)). Notice that the upper most sub-band, 
laying closest to the Fermi energy, is well separated in 
energy from the rest and has a bandwidth of about W ~ 
45meV. The Coulomb on-site energy U is evaluated to 
be on the order of O.leV [10 . Therefore, the low energy 
physics can be described by a single band model with the 
value of U/W ^ 2.2 well over the critical value for the 
Mott transition(/7c/W ^ 1.3) [16 . 

Finally, we see a large hysteresis in the resistivity 
measurements as a function of temperature as shown in 
Fig.|2|g). This is characteristic of a first order phase tran- 
sition into the CCDW phase. Disorder changes somewhat 
the transition temperature and the hysteresis. Collec- 
tively, all the experimental data suggests that the com- 
mensurate inter-locking of the "star-of-David" clusters in 
the CCDW phase is not affected by Cu intercalation. 

We next investigate the effect of disorder on the low 
energy spectral properties. In Fig[3]^a-b) we show a com- 
parison of APRES data of a clean and a disordered sam- 
ple, at the same crystal orientation along the F — di- 
rection measured with 22 eV photon energy. There is a 



significant increase in the spectral intensity in the disor- 
dered sample in the [—0.1 eV, eV] energy range which 
appears to be centered at the F point. Looking at the 
EDCs at the F point in Fig [sjc) one can see the trans- 
fer of spectral weight towards lower binding energy in 
the disordered sample, with a line shape that is slightly 
skewed towards E=0. The symmetrized EDC (Fi^Sjd)) 
confirms the presence of a "soft" gap, where the spectral 
weight at zero energy is about 20 % of the intensity at the 
peak and the broad peaks are shifted inwards, leaving a 
smaller Mott gap. The filling up of Mott gap is accompa- 
nied by a transition from insulating to metallic behavior 
in resistivity (Fig. |2|g)). However, it is a very unusual 
disorder driven metal with no Fermi liquid like spectral 
peak at cj = 0. In Fig. [sj^e) we plot the EDCs at the F 
point for 13 different clean and disordered crystals. The 
results are highly reproducible, although the data shown 
here were measured in different ARPES setups and the 
crystals are from several different growth batches. 

The temperature dependence of the symmetrized EDC 
at F is shown in Fig. [sj^f) for a clean and disordered 
sample. The difference in line shape and in the shape 
of the gap is most prominent at low temperatures, and 
the EDCs become increasingly similar as temperature is 
raised. 

To gain better insight into our experimental results, we 
investigate the effect of disorder on a Mott insulator using 
a half-filled single band Hubbard model on a triangular 
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FIG. 3: The electronic structure of a disordered Mott-insulator: (a) ARPES data from a lT-TaS2 crystal taken with 22eV 
photons along the T — K direction at T< Tccdw, reveal a Mott gap in the spectrum as indicated by the lack of intensity 
within about 80 meV of the Fermi energy, (b) ARPES data from a disordered samples measured under the same conditions; 
significant amount of spectral weight is now present close to the Fermi energy, inside the energy gap. (c) Comparison of 
the EDC at the F point taken from the scans in (a) and (b). The peak in the disordered lT-TaS2 is broadened and shifted 
towards lower binding energy and there is substantial amount of intensity at zero energy, (d) The same EDCs shown in (c) 
after symmetrization, clearly reveals the closing of the gap in the disordered sample, the intensity at the Fermi-level is about 
20 % of the peak- intensity, (e) Reproducibility of the results: EDCs from 13 different samples (7 clean and 6 disordered) 
taken in different ARPES systems with different photon energies and polarizations, all reveal the same qualitative behavior, (f) 
Symmetrized EDCs at the F point as function of the temperature. The samples were first cooled to base temperature and 
measured while warming to room temperature. While at low temperatures the spectral weight at zero energy is much larger 
in the disordered samples, above the CCDW transition temperature the spectral weight at zero in the clean and disordered 
sample is the same. 



lattice. It describes the CCDW phase of IT-TaSs- 

(*i>,cr {{ij)) 

Here t is the nearest-neighbor hopping between unit cells, 
t' the next nearest hopping, /i, the chemical potential 
and [/, the on-site interaction strength. c\^{cia-) cre- 
ates (annihilates) an electron of spin <j at site and 
^icT = ^Icr^io- is the numbcr operator for spin a at site 
i. The onsite disorder potential e^, with a flat distri- 
bution in the range [— y/2,V/2], captures the effect of 
lattice distortions due to Cu intercalation. In our analy- 
sis, addition of Cu does not change the filling. We solve 



Eq.[T] within inhomogeneous mean- field theory where in- 
dependent mean-field parameters at each site are deter- 
mined self-consistently (details in Supplementary Mate- 
rial). For a given disorder realization, we obtain the 
single particle eigenstates {|^ck)} and eigenvalues {coc} 
which can be used to calculate the single particle Green's 
function 

G(ri,r2,c.) = yi^:^^^^^K^ (2) 

a 

Upon disorder averaging, and then Fourier transform- 
ing, we finally get the spectral function A(k,uj) = 
— (l/7r)Im G(k, cj) which we compare with experimental 
results below. 

EDC of spectral function: At half- filling in the clean 
limit with f = 0, mean field solutions gives an insulator 
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FIG. 4: (a) Theoretical density of states for different disorder 
strengths V and fixed interaction JJ — 6t. (b) EDO of spectral 
function k — kp for U — and different values of disorder 
strength V . For weak disorder (!/<[/), there is a well formed 
gap which closes as disorder is increased {V ^ U). Inset: 
Spectral weight at cj = as a function of disorder strength. 
For V there is spectral weight at cj = 0. We used 18x18 
size lattices for the calculations averaged over 100 disorder 
realizations, (c) and (d) shows probability density of a; = 
eigenstate as a function of position for U — and U — 
respectively with V — 7t. The eigenstate is extended over the 
entire lattice for U — 6t. The scale is cut off at the maximum 
value of 1^(01^ in each case. 



for U > Uc — 5.27t [17 along with magnetic ordering 
of the spins on the triangular lattice in a commensurate 
120*^ pattern. For the rest of the calculation, we choose 
U = 6t in order to remain in the insulating phase. In- 
clusion of t' = 0.3t does not affect the magnetic ground 
state, but is important for producing a realistic band 
structure for comparison with experiments. Small change 
in filling {5n < 2%) also does not change the mean-field 
solution. 

Previous studies have observed a suppression in the to- 
tal density of states around u = persisting up to large 
values of disorder jT8H2Q] . Here we obtain information 
about the spectral function in Fig. ^h). The EDO at 
the Fermi momentum for different amounts of disorder 
can be compared directly with the experimental results 
shown in Fig. |3jc-d). Our theoretical results agree quali- 
tatively with our experimental results on several aspects. 
First of all, with increasing disorder the Mott gap closes 
and finite spectral weight develops at = (see inset of 
Fig. [4]^b)). Second, the shape of the gap changes from a 
hard 'U' shaped gap in the clean limit to a 'V shaped 
soft gap in the disordered case. Finally, the position of 
the broad peaks corresponding to the upper and lower 
Hubbard bands move closer to the chemical potential 
with increasing disorder in agreement with experimental 



data. It should be pointed out that the spectral function 
from theory is asymmetric around u = and is due to 
the underlying asymmetry of the density of states (DOS) 
shown in Fig.Qa). The experimental data has been sym- 
metrized around = to remove the effect of the Fermi 
factor for small energies (cj <C A where A is the gap 
between the broad peaks). Comparison between theory 
and experiment should be made only for small cu. 



In our calculation, the additional low frequency weight 
upon disordering the system arises primarily from k- 
states close to the underlying Fermi surface of the tight- 
binding model. In the experimental data in Fig. |3|b), 
on the other hand, the additional weight appears from 
states close to the F point. The discrepancy could be 
due to experimental resolution or possibly due to mag- 
netic ordering and the resultant zone folding that can 
also push the weight close to the T point. 



Mid- gap states: The crucial question to ask is about 
the nature of states filling up the Mott gap when dis- 
order is introduced. The scaling theory of Anderson 
localization for non interacting electrons predicts that 
these states should be localized, however, the interplay of 
strong interaction and disorder can lead to a novel metal- 
lic phase [9 . We confirm here that the eigenstates close 
to the chemical potential are indeed extended on a trian- 
gular lattice for V ^ U (See Fig. |4jc-d)). This is further 
confirmed by the observation of a drastic drop in resis- 
tivity of Cu intercalated lT-TaS2 in the CCDW phase 
in comparison to the clean samples (Fig. |2|g)), thereby 
providing evidence for the emergence of this novel metal. 



In conclusion, the remarkable "star-of-David" recon- 
struction and the strong correlations on the Ta 5d or- 
bitals conspire to make lT-TaS2 a realization of a single 
orbital Mott insulator on a triangular lattice. We have 
found that lT-Cua;TaS2 with small values of x (less than 
1%) does not destroy the underlying CDW state showing 
crucially that Cu intercalation disorders the system but 
does not add carriers. Using ARPES we have shown that 
disorder destroys the Mott gap and introduces electronic 
states at zero energy. Based on the comparison between 
the theoretical results and the ARPES and transport ex- 
periments, we conclude that lT-Cua;TaS2 is the first ex- 
ample of a novel metallic phase in a disordered 2D Mott 
insulator. This conducting state is a non-Fermi liquid 
metal with a pseudogap in the spectral function. Our in- 
vestigations open up the possibility of several future ex- 
plorations of the nature of magnetism, its evolution with 
disorder and its relation to the metal-insulator transi- 
tions. 
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METHODS 
Sample Preparation 

Single crystals of lT-TaS2 were prepared using the va- 
por transport method as described in Ref. [22j. For the 
addition of disorder, we added 10 % Cu into the stoichio- 
metric mix during the crystal growth stage. However, the 
copper content of the single crystals that emerge from the 
growth is well below 10 %. In fact, wave dispersive elec- 
tron spectroscopy places the Cu content of the sample 
at 0.75% ± 0.05%. The Cu atoms are expected to re- 
side as intercalation between the TaS2 lavers[2T]. These 
samples, which we shall refer to as disordered lT-TaS2, 
display a lower NCCDW-CCDW transition temperature 
of around 150 K upon cooling (Fig.|2|g)), and the nega- 
tive slope of R{T) is also considerably reduced compared 
to the clean sample, suggesting enhanced conductivity. 

ARPES 

We have performed extensive angle-resolved photoe- 
mission spectroscopy (ARPES) measurements, on the 
PGM-beam line at the Synchrotron Radiation Center 
with 72 eV and 22 eV photon energy as well as at the 
ARPES machine at the Technion, using the Hel 22 eV 
spectral line from discharge lamp as the light source. In 
both cases data was collected using a scienta R4000 elec- 
tron analyzer. Samples were cleaved in-situ at a pressure 
better than 5 x 10~^^ torr. Energy and angular resolution 
were better than 10 meV and 0.3°. 
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SUPPLEMENTARY MATERIAL 
Disorder effects in the NCCDW 

We verify that the new states reported in the paper 
are linked to the underlying Mott state by searching for 
any sign of increased spectral weight near the Fermi level 
in the normal metal phase that emerges from the non- 
commensurate charge density wave (NCCDW) phase. In 
Fig. [l] we compare the EDCs at the F point and at k^? of 
clean and a disordered sample. The samples were cleaved 
and measured at 300 K and the data normalized by the 
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FIG. 1: EDCs of a clean and disordered sample taken with 22 eV photon energy at room temperature at F in (a) and at kp 
in (b). 



intensity at —0.8 eV. The EDO's are identical and this 
remains the case throughout the Brillouin zone. We find 
no trace of the new states when the disordered sample 
is not in the Mott state. In addition, we measured very 
carefully the position of the top of the Sulfur band with 
respect to the Fermi-level. The top of Sulfur band at 
about 1.2eV below the Fermi- level is relatively sharp and 
can be used to measure the chemical-potential. We found 
no measurable effect of Cu intercalation on the chemical- 
potential. We conclude that if indeed the Cu intercala- 
tion induces new states, these states cannot account for 
the additional spectral intensity within the Mott gap in 
the commensurate CDW (CCDW) phase. 



Inhomogeneous mean field theory 

The commensurate CDW (CCDW) phase of IT-TaSs 
is widely regarded as a realization of half-filled 2D Hub- 
bard model on a triangular lattice [1-3 . Of the 13 Ta 
5d electrons in the "star-of-David" unit cell 12 occupy 
states below the Fermi energy, while the 13th electron 
occupies a half-filled band crossing the Fermi energy. In 
the case of Cu intercalated samples, we need to include 
the effect of disorder also, which we model as on-site po- 
tential disorder. The model Hamiltonian is (Eq. (1) in 
main text) 

(*i>,cr {{ij)) 

+ ^ (ei - /i) 4^Q^ + ^ XI (1) 



TEM 



Transmission Electron Microscopy (TEM) probes the 
periodic structure in the material and is ideal for de- 
tecting the presence of the electronic CDW modulation 
in the sample, and identifying the different phases of 
TaS2. The CCDW-NCCDW transition can be seen in 
the clean sample in Fig. |2] (top panels), above the tran- 
sition the brightest spots correspond to the diffraction 
from the undistorted hexagonal Ta sheets, and the larger 
dim hexagons are associated with the formation of the 
" star-of-david" unit structure. Below the transition the 
stars interlock and tile the entire Ta sheet, yielding the 
new diffraction peaks below the transition. We observed 
exactly the same behavior in the disordered samples, as 
shown Fig. |2] (bottom panels). The addition of Cu did not 
suppress the CCDW phase. The electron diffraction im- 
ages were acquired using a FEI Tecnai T20 Transmission 
Electron Microscope (TEM), using a Gatan double-tilt 
cryo-stage (liquid N2). The samples were cooled down to 
lOOK and then warmed up back to room temperature. 



Here t is the nearest-neighbor hopping between unit cells, 
t' the next nearest hopping, /i, the chemical potential 
and [/, the on-site interaction strength. c\^{cicr) cre- 
ates (annihilates) an electron of spin cr at site z, and 
^icr — ^ia^icr IS the uumbcr operator for spin a at site 
i. The onsite disorder potential e^, with a fiat distribu- 
tion in the range [— V/2,]//2], captures the effect of lat- 
tice distortions due to Cu intercalation. In our analysis, 
addition of Cu does not change the filling. 

We solve Eq. [T] within inhomogeneous mean-field the- 
ory. The interaction term is approximated as 



-{st)s--{s-)si^{st){s-) 



(2) 



where the operator S = ^> c^Tcj^ct'Ccj' and r are the 
Pauli spin matrices. Due to the presence of disorder, the 
density and spin fields are site-dependent and the 3A^ pa- 
rameters ({(rii^), = {S^Y}, i = 1, . . . , TV and cr =t 
or ^), where N is the total number of sites, are obtained 
by a self-consistent solution of the coupled inhomoge- 
neous Hartree-Fock equations. Once self-consistency is 
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FIG. 2: TEM diffraction images of a clean and disordered samples in the CCDW phase (left panels) and in the NCCDW phase 
(right panels) 





FIG. 3: (a)&(c): Momentum distribution cut (MDC) of disorder averaged spectral function at cj = plotted in the magnetic 
Brillouin zone (dotted hexagon) for = and t' = 0.3t respectively. (b)&(d): The same MDC plot as in (a) and (c) but in 
the non-magnetic Brillouin zone. V=7t and U=6t 



reached, we obtain the single particle eigenstates and 
eigenvalues, and use them to calculate the spectral func- 
tion. For a given disorder realization, the single particle 
Green's function is given by (Eq. (2) in main text) 



G(ri,r2,cj) = ^ 



{yi\^Iji){^Iji\v2) 
CO ^if] - ei 



(3) 



where is the mean- field eigenstate with eigenvalue e^ 



Upon disorder averaging, we obtain the translationally 
invariant Green's function, G(ri,r2,co') G{ri — r2,co') 
which can be Fourier transformed to get the Green's func- 
tion in momentum space, G(k, a;). The spectral function 
can finally be calculated as 

A(k,cj) = --ImG(k,cj) (4) 

TT 

Magnetic order: Inhomogeneous mean-field theory 
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predicts magnetic order with the commensurate 120^ pat- 
tern even in the presence of disorder. While there is 
no direct experimental evidence of magnetic order, the 
APRES data shown in Fig. (3) of main text show sig- 
nificant uj = spectral weight at the F-point. We can 
reproduce the spectral weight at F-point from theory 
only if we include the zone folding effect of the magnetic 
lattice in /c-space. Fig. jsjb) and (d) show the momen- 
tum distribution cut (MDC) of the spectral weight at 
= for two different values of t' (next-nearest hopping 
amplitude) plotted in the non-magnetic Brillouin zone. 
There is negligible spectral at the F-point. In compari- 
son, Fig. [sj^a) and (c) which are plotted in the magnetic 
Brillouin zone (dotted hexagon) show significant spectral 
weight at the F-point. The observed spectral weight at 
the F point suggests enhanced magnetic correlations at 
the commensurate ordering wave vector. We predict that 



experiments performed at lower temperatures should ob- 
serve long range magnetic order. 
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